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1. INTRODUCTION 

In recent years, ever increasing energy demand has increase the stress on power system which 
eventually would cause a power system to operate near its stability limit [1, 2]. In order to mitigate this 
situation, several scheme can be implemented to improve the stability of a power system such as installation 
of power generation unit [3], reconfiguration of power networks [4] as well as installation of compensation 
devices. 

Flexible AC Transmission System (FACTS) devices can be installed to provide compensation to 
power system. Static Var Compensator (SVC) is a type of shunt-connected FACTS devices which is capable 
of injecting and drawing reactive power from or to a power system. Due to this flexibility, SVC is preferred 
over capacitor bank in providing reactive power support to a power system [5]. However, optimal installation 
of SVCs are important for an effective compensation. To achieve this optimality, optimization techniques are 
implemented to determine the optimal location and sizing of SVCs to be installed in a power system. 
Various studies have been conducted to optimally allocate FACTS devices in power systems to solve various 
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issues related to power system operations. Authors in [6] has determined the optimal allocation of SVC on 
IEEE 30-Bus using Genetic Algorithm (GA). Similar optimization technique is also implemented by [7] to 
determine the optimal allocation of SVC on a 13-Bus power system with the aim of improving voltage 
profile, minimize power loss and minimize investment cost. In [8], Raj and Bhattacharyya has implemented 
Whale Optimization Algorithm (WOA) to optimally allocate SVC and TCSC in IEEE 30-Bus and IEEE 57- 
Bus systems to solve reactive power planning problem. An attempt to minimize active power loss and 
improve voltage profile of 57-bus test system via optimal allocation of shunt capacitor, SVC and STATCOM 
has been conducted by Sirjani et al. [9] by using Novel Global Harmony Search (NGHS) algorithm. 
Bhattacharyya and Gupta [10] has solve reactive power optimization problem using fuzzy-based evolutionary 
algorithms in which that the location of the FACTS devices are determined using Fuzzy Logic while the 
evolutionary algorithms will determine the optimal settings for the power system. The evolutionary 
algorithms used in [10] are GA and Differential Evolution (DE). A modified version of Gravitational Search 
Algorithm (GSA) known as Cumulative Gravitational Search Algorithm (CGSA) has been proposed to solve 
optimal STATCOM allocation problem in power system [11]. 

Although metaheuristics-based optimization techniques are able to solve optimal FACTS device 
allocation problem, problems have been reported related to the performance of the optimization technique 
such as dependency to parameter of optimization technique, tendency to stuck in local optima, 
slow convergence time and high computational burden, which eventually prevent the optimization technique 
from producing high quality solutions [5, 12]. Hence, a higher-performance optimization technique is needed 
so that the FACTS device can be optimally installed in a power system which would provide effective 
compensation. Symbiotic Organisms Search (SOS) optimization technique was first developed by Cheng & 
Prayogo [13] in 2014 which was inspired by the symbiosis relationship of organisms in the ecosystem. 
The performance of SOS has been acknowledged by several studies conducted. In previous studies, SOS has 
been implemented to solve optimal SVC allocation for voltage profile improvement [5], economic dispatch 
[14, 15], optimal power flow [16] as well as civil engineering design [17]. Further modifications are 
proposed on SOS to improve its performance. Tejani et al. [18] has introduced adaptive element on the 
benefit factor of the SOS. Abdullahi et al. [19] has substituted the uniformly-distributed random number on 
the mutualism and commensalism updating expression by using random number generated by using chaos 
sequence. Saha et al. [20] has introduced chaotic local search (CLS) algorithm before the convergence test 
which would attempt to find a better solution around the current best solution. This methodology has been 
implemented and validated by Mohamad Zamani et al. [21]. 

This paper presents the implementation of Chaotic Immune Symbiotic Organisms Search (CISOS) 
optimization technique to solve optimal SVC allocation for voltage security control. In this paper, the 
location and sizing of the SVC will be determined by CISOS. The optimization problem will be 
demonstrated on IEEE 26-Bus Reliability Test System (RTS). Comparative studies are also conducted with 
SOS, Particle Swarm Optimization (PSO) and Evolutionary Programming (EP) highlights the superiority of 
CISOS in solving the optimal allocation problem in terms of solution quality produced by the optimization 
techniques. 


2. RESEARCH METHODOLOGY 

In order to solve the optimal allocation of SVC problem for voltage security control, CISOS 
technique is applied to determine the optimal location and sizing of SVC to be installed in the test system to 
improve the voltage security of the system. 


2.1. Problem Formulation 

The first objective of the study is to improve the voltage profile of the system through the 
minimization of total voltage deviation index (TVDJ). The fitness function for the first optimization objective 
fi as well as mathematical representation of TVDIJ are expressed as in (1). 


2 
Vre fini 
(iS hVDray (=) (1) 


Vref,ji 


where V,.¢; is the reference value of i” bus voltage, V; is the i” bus voltage and N, is the number of buses in 
the test system. In [5], TVDI is calculated only for load buses. For this paper, the calculation of TVDIJ will 
include all buses in the system. The voltage reference value for load buses are set to 1.00 p.u. while the 
reference value for slack and P-V buses are set with their voltage setpoint value as in the power flow solution 
algorithm [21]. 
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The second objective of the study is to improve the voltage stability of the system. The voltage 
stability of the system is evaluated by using Fast Voltage Stability Index (FVS/J) developed by Musirin et al. 
[22]. Voltage stability of the system is represented by the highest FVSI value in the entire test system. 
FVSI is mathematically expressed as in (2) while the fitness function of the second optimization objective f> 
is expressed as in (2). 


422.40 
FVSlqp = Tae (2) 
fo = FVSI max = Max(FVSIq,) (3) 


where FVSI ng, is the highest FVSI value from the entire test system and FVSI,, is the FVSI value for the 
transmission line connecting a” bus and b” bus in the test system. The sending bus and receiving bus of the 
transmission lines are represented by a” bus and b” bus respectively. Voltage value of the sending bus is 
represented by V, and the reactive power flowing into the receiving bus is represented by Q,. The reactance 
and impedance of the transmission line connecting a” bus to b” bus are represented by X,, and Za, 
respectively. 

Normally, optimization techniques are used to solve single-objective optimization problem. 
However, this study requires the proposed optimization technique to solve the optimization problem 
considering both objectives. In order to achieve the goal, multi-objective optimization mode should be 
implemented. In this paper, weighted sum method is applied to solve multi-objective optimization problem 
by combining all fitness functions and weighted by using weighing factor [23]. Hence, the objective function 
of the optimization is expressed as: 


0. F.= minimize(f) (4) 


f=, x fi) + 2 x fr) (5) 


where f is the overall fitness value which is the sum of the weighted fitness values and O.F. is the objective 
function. The weighing factor for the first and second fitness values are represented by w; and w>. 


2.2. Constraints 

In weighted sum method, sum of the weighing factors used in (5) should be equal to 1 [24]. 
This constraint can be expressed as in (6) where w; is the weighing factor of i" fitness function and Nouj is the 
number of optimization objectives. 


Dew = 1 (6) 

It is possible for an SVC to inject or draw reactive to or from a bus in a power system. In this paper, 
the SVC is modelled as a negative reactive power load. Positive SVC sizing value Q,,, indicates the reactive 
power from the SVC into the bus. On the other hand, negative value of Q,,. implies the value of reactive 
power drawn from the bus into the SVC. The sizing of SVC is limited by its minimum and maximum 
capacity of the SVC indicated by Q™” and Q™4* respectively. This constraint can be expressed as: 


Ver S0.9 5 On (7) 


In this paper, SVC is proposed to be installed at load buses since slack bus and P-V bus have the 
capability to regulate reactive power. Load bus is defined as a bus with a load connected to it. A bus with no 
load or power generation unit connected to it is known as connecting bus, and SVC is not proposed to be 
installed there. The constraint of location of SVC installation can be expressed as: 


1<LOC,,. < Ny (8) 
where N;, is the number of load buses in the system. 
2.3. Chaotic Immune Symbiotic Organisms Search (CISOS) for Optimal SVC Installation 

To solve the optimal SVC allocation problem, the authors have proposed the implementation of 


CISOS optimization technique to determine the optimal location and sizing of SVC for voltage security 
control in power system. CISOS is first developed by Mohamad Zamani et al. [25] which is inspired by the 
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inclusion of diversification and improved exploration capability through the introduction of element of 
cloning, elimination and Chaotic Local Search (CLS) to the original Symbiotic Organisms Search (SOS) 
optimization technique developed by Cheng and Prayogo [13]. Detailed information about CISOS can be 
referred in [25], and the process of CISOS in solving optimal SVC allocation problem for voltage security 
control can be summarized as follows: 


2.3.1. Step 1: Initialization 

Firstly, the optimization parameters are set which are number of initial organisms in the ecosystem 
Nj, Number of clones Neione, iteration counter iter, maximum number of iteration max_iter, and CLS control 
parameter p. Then, a set of possible solution known as the organism is generated by using random number. 
A set of organism is made up of 2 parts, namely SVC sizing and the location of SVC. The SVC sizing is 
generated using random number bounded with the limits specified as in (7) while the SVC location is 
generated using integer random number bounded by limits specified as in (8). 


iter =1 (9) 
Osveaay = rand(Qhn, One" (10) 
LOC yea) = randInt(1, Nip) (11) 


Qsve(1,1)) oy Qsve(1,d)» Fey Qsvc(4,Neve)? LOC sy(1,1) aie LOC sye(1,4)) ek) LOC sy6(4,Neye) 
Xinit = Qsve(i,1)» sey Qsve(i,d)s Ty QsvcliNeve)? LOC syc(i,1)) a LOC syc(i,a)s 7) LOC} S46 (i,Neye) (12) 
Qsve(Wio)? my Qsve(Wiod)? aa Qsve(WioNeve)? LOC sy¢(Wigs)? ee LOC syc(Nig.d)) ar LOC sye(Nig Neve) 


while Q,,.;,a) is the sizing of d" unit SVC of i” initial organism, LOC,y¢ii,2) is the location of d unit SVC of i” 
initial organism, N;, is the number of initial organisms in the ecosystem and N,,, is the number of SVC to be 
installed. In (10), rand() is the function to generate random number and randInt() in (11), is the function to 
generate integer random number. 

After the ecosystem has been occupied by the initial organisms, then the fitness value of the initial 
organisms is computed. Then, an organism with the best fitness is chosen to be initialized as the best 
organism. Lastly, the chaotic search space radius is computed by using (13), which consists of (14) and (15). 


= {rye TG.) TAG Neve) T2090 T24,a)0 T2Ngve)} (13) 


Orretd) ~ eve) 
SVC SVC 
"1aa)-=— 3 (14) 


max min 
LOC sy¢(d) LOC Sy 6(d) 


T2(,4) = 5 (15) 


2.3.2. Step 2: Cloning 

Size of the ecosystem is increased by cloning the initial organisms. The initial organisms will be 
duplicated by a factor of Neione. After the cloning process, the organisms X will be subjected to the next 
optimization step. The act of cloning to produce X can be expressed as: 


(2) 
x=| Xinic (16) 


x Netone) 


init 


After cloning, the number of organisms in the ecosystem has increased by a multiple of Netone, and 
the total number of organisms in the ecosystem after cloning N,7 can be computed as: 


Nor = Nio x Netone (17) 
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After the initial organisms have been duplicated, then the organism counter i is initialized to 
value of 1. 


i=1 (18) 


2.3.3. Step 3: Symbiosis Relationship Phases 

The organisms are subjected to symbiosis relationship phases consisting of mutualism, 
commensalism and parasitism. Firstly, an organism known as the i” organism X; is chosen from the 
ecosystem. Then, another organism known as j” organism X; is randomly chosen from the same ecosystem. 
When choosing X,, ” organism should not be the same organism as the i” organism. Then, both organisms 
are subjected to mutualism phase, which then produce new organisms represented by Xinew and Xjnew. Xinew 
and Xj new are produced by 


Pa Xj+Xj 


MV = 5 (19) 
BF1 = randint(1,2) (20) 
BF2 = randint(1,2) (21) 
Xinew = X; + rand(0,1) X (Xpese — MV x BF1) (22) 
Xj new = Xj + rand(0,1) X (Xpese — MV x BF2) (23) 


where MV is the mutual vector computed by (19), X45; is the best organism in the ecosystem, BF/ and BF2 
are benefit factors, rand(0,1) is a random number ranged between 0 to 1, and randInt(/,2) is an integer 
random number between | to 2. It should be noted that location of SVC in Xj ney and Xjney produced from the 
updating process should be rounded off to the nearest integer. 

Then, the fitness values of Xj,ne, and Xj, are computed and compared with the fitness values of X; 
and X; respectively. If the fitness value of Xjne, is lower than X;, then X;,,., will replace X; in the ecosystem. 
The same case also applies to Xj, and Xj. 

Next, the organisms are subjected to commensalism phase. Maintaining the same X;, another random 
organism is chosen from the ecosystem to be X; in which that X; is not the same organism as X;. Then, a new 
organism X;,ne, 1s produced by subjecting X; and X; to commensal relationship which is expressed as 


Xinew =X; + rand(—1,1) x (Xpest — Xj) (24) 


where rand(-1,1) is a random number between -1 to 1. The location of SVC in X;,2, produced by (24) should 
be rounded off to the nearest integer. The fitness value of X;,,2, 18 computed and compared with the fitness 
value of X;. Should X;,2, possess a lower fitness value compared to X;, then X;,,., will replace X; in the 
ecosystem. 

Final phase in symbiosis relationship phases is parasitism phase. The same X; is maintained, 
an organism X; is randomly chosen from the ecosystem. The same rule for choosing X; in mutualism and 
commensalism also applies to choosing X; in parasitism. Then a parasite PV is produced by duplicating X;. 
Then random dimension of PV is modified by using random number as in (10) for SVC sizing and (11) for 
SVC location. Then, fitness value of PV is computed and compared to fitness value of X;. Should the fitness 
value of PV be lower than X;, then PV will replace X; in the ecosystem. 


2.3.4. Step 4: Best Organism Identification 

A candidate for best organism is chosen among the organisms in the ecosystem in which that the 
organism possesses the lowest fitness value. Then the fitness value of the candidate is compared to the fitness 
value of current best organism. If the best organism candidate has a lower fitness value compared to the 
fitness value of current best organism, then the best organism candidate will replace the current best 
organism as the best organism. Otherwise, the current best organism is retained. 


2.3.5. Step 5: Symbiosis Relationship Termination Test 
The optimization algorithm will check either all the organisms has undergone the symbiosis 
relationship phases. If the organism counter i has reached total number of organisms after cloning N,7, 
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the optimization process will continue to step 6. Otherwise, the organism counter is increased by 1, and the 
optimization process continues to step 3. 


2.3.6. Step 6: Elimination 

The size ecosystem is reduced from N,7 organisms to Nj. To reduce the ecosystem size, CISOS will 
choose Nj, fittest organisms among the entire organisms in the ecosystem. This way, only the fitter organisms 
will be carried forward to next iteration. 


2.3.7. Step 7: Chaotic Local Search (CLS) 

A new organism is generated around the best organism using CLS. Firstly, CLS iteration counter ch 
and chaotic variable cv are initialized. The number of dimensions of the chaotic variable should be equal to 
the number of dimensions of Xpes;. 


ch =0 (25) 
CV(ch,c) = rand (0,1) (26) 


where c is the dimension of cv. 
Then, chaotic variable is updated by using Piecewise Linear Chaotic Map (PLCM). 


CV(ch,c) 


CVn ey H1OD) 


CV(ch+1,c) = (1-cv¢en,c)) ” 
aap eM che) © (PD) 


where p is the control parameter of CLS. 

After the chaotic variable is updated, a new organism is produced by wrapping the chaotic variable 
around the best organism X;,s, using (28). It should be noted that the location of SVC produced by (28) 
should be rounded off to the nearest integer. 


X cls(ch+1,c) = X pest(1,c) an Fea 2606s ~e 1) (28) 


The fitness value of the produced organism is computed and compared with the fitness value of the 
best organism. If the fitness value of the new organism is higher than the fitness value of Xj.5,, then the CLS 
iteration counter is increased by | and the CLS algorithm proceeds to chaotic variable update. 


ch=ch+1 (29) 


In the case of new organism having lower fitness value compared to X45, then the new organism 
will replace X;,.s,, CLS algorithm is terminated and the chaotic search space radius is reduced by using (30). 
In (30), rand(0,1) is arandom number ranged from 0 to 1. 


T(1,c) = (1c) + rand(0,1) (30) 


2.3.8. Step 8: Convergence Test 

In this step, CISOS will determine whether the algorithm should continue or stop. If the value of iter 
has reached the maximum CISOS iteration number max_iter, then the optimization process will stop, 
and X;,.5, will be the optimal solution produced by CISOS. Otherwise, value of iter will be incremented by 1, 
and the optimization process continues to step 2. 


3. RESULTS AND DISCUSSION 

In this paper, CISOS is implemented to solve optimal SVC allocation for voltage security control 
problem in power system. The test system used in the study is IEEE 26-Bus Reliability Test System (RTS). 
IEEE 26-Bus RTS consists of 17 load buses and 6 power generation units. In the optimization process, 
3 units of SVC will be installed in which that the location and sizing of the SVC will be determined using 
CISOS. The parameters of the optimization process are listed in Table 1. 
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Table 1. Optimization parameter of CISOS 








Optimization parameter Parameter value 
Number of generated organisms 20 
Number of clones 10 
Maximum iteration of CISOS algorithm 200 
CLS control parameter 0.4 
Maximum CLS iteration 100 





The optimization process will run for 20 times for each case study to observe the variation of 
optimal solutions produced by CISOS. The results which would be analysed is the fitness value of the 
optimal solutions which are made up of TVDI and FVSInq, values. The effectiveness of CISOS in solving the 
allocation problem is demonstrated through comparative studies with respect to Particle Swarm Optimization 
(PSO) and Evolutionary Programming (EP) techniques. 


3.1. Base Case Condition 

During this condition, the test system experience no change from its nominal condition. 
Optimization techniques are implemented to solve the optimal SVC allocation problem on IEEE 26-Bus 
RTS. The best optimal solutions produced by CISOS, PSO and EP are tabulated as in Table 2. 


Table 2. Results of Voltage Security Improvement Using SVC During Base Case Condition 
Best Optimal Solutions 








Fechuidue SVC Sizing (MVAr) SVC Location (Bus No.) _ TVDI FVSInax Fitness Value 
Pre-optimized = - 0.00483 0.35376 0.17929 
CISOS 60.2433, 40.9218, 37.5464 9,12, 11 0.00463 0.33607 0.17035 
PSO 100.0000, -76.6738, 100.0000 6,9, 22 0.00162 0.34324 0.17243 
EP -0.8570, 71.2966, -86.8921 21, 20, 18 0.00292 0.34294 0.17293 





From the results obtained, it can be observed that the optimization techniques are able to solve the 
optimal SVC allocation problem. All techniques have produced optimal solutions with lower fitness values 
compared to the pre-optimized fitness value. The results also revealed that CISOS has managed to record the 
best optimal solution which having the lowest fitness value compared to the optimal solutions produced by 
PSO and EP. The optimal solution produced by CISOS has improved the voltage profile of the system 
indicated by the reduction of TVDI value from 0.00483 to 0.00463. However, the optimal solution produced 
by CISOS provides maximum improvement to the voltage stability of the system which is indicated by 5% 
reduction of FVSIingx value from 0.35376 down to 0.33607. The optimal solution produced by EP provides 
better voltage profile improvement, but it only manages to produce lesser improvement on voltage stability of 
the system. Although the optimal solution produced by PSO provides maximum improvement on the voltage 
profile of the system, the improvement on the voltage stability is minimal compared to CISOS and EP. 


4. CONCLUSION 

This paper has presented the implementation of Chaotic Immune Symbiotic Organisms Search 
(CISOS) optimization technique to solve optimal SVC allocation problem on IEEE 26-Bus RTS for voltage 
security control. It is desirable for a power system to have a good voltage security since a good voltage 
security suggests a power system to be healthy and able to maintain stable operation. Implementation of 
optimization techniques can help to achieve this while mitigating overcompensation and under-compensation 
issues. While the fitness value has been reduced through the implementation of optimal solutions produced, 
the superiority of performance of CISOS compared to PSO and EP has been proved in terms of the quality of 
the optimal solution produced. 
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